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I .  INTRODUCTION 


The  interaction  of  a  turbulent  wind  trlth  a  dust-laden  surface,  and  the 
subsequent  lofting  of  dust  particles,  is  of  interest  for  a  variety  of 
applications . 


If  the  surface  particles  are  relatively  heavy  and  the  wind  speed  is  not 
large,  the  individual  particles  ejected  from  the  ground  rise  a  certain 
distance,  travel  with  the  wind,  and  then  descend  back  to  the  ground.  This 

process,  referred  to  as  "saltation,*'  leads  to  the  generation  of  desert 

12  2 
dunes.  *  The  saltation  region  occurs  in  the  range 


-2  K 

0(10  )  <  ^  5  rrr-A 
F  og  d 


<  0(1) 


(1) 


where  Fq  is  the  aerodynamic  drag  force  on  a  particle  resulting  from  turbulent 
fluid  motion  and  F.  is  the  gravitational  force.  Other  symbols  are  defined  in 

O 

the  Appendix.  The  lower  limit  is  related  to  the  minimum  force  required  to 
loft  a  particle,  and  the  upper  limit  corresponds  to  the  limit  beyond  which  the 
particle  goes  into  suspension.  For  sand  particles  in  standard  air,  Eq.  (1) 
becomes  approximately  O(IO^)  <  u_^/d  <  0(10  ),  where  u^  (cm/sec)  is  the  local 
friction  velocity^  (l.e.,  the  characteristic  turbulent  eddy  velocity  near  a 
wall)  and  d(cm)  is  particle  diameter.  For  d  ••  10  ^  cm,  the  friction  velocity 
is  limited  to  values  between  10  cm/sec  and  100  cm/sec  in  order  for  saltation 
to  occur. 


With  increase  in  wind  velocity,  the  dust  particles  enter  suspension. 

This  process  is  of  interest  in  the  generation  of  dust  clouds  and  in  the 
evaluation  of  nuclear  weapons  effects.  Relatively  few  experimental  studies 
have  been  made  of  dust  ingestion  into  a  high-speed  turbulent  boundary  layer. 
Hartenbaum^  has  Investigated  air  flow  over  a  semi-infinite  sand  bed  with  free- 
stream  velocities  in  the  range  3  x  10^  to  10^  cm/sec  and  a  mean  particle 
diameter  d  -  0.025  cm.  Ausherman^  has  investigated  the  turbulent  boundary 
layer  behind  a  shock  moving  over  a  sand  bed  with  a  mean  particle  diameter 
d  ■  0.009  cm.  Flow  velocities  relative  to  the  wall  were  in  the  range  2  x  10^ 
to  5  X  10^  cm/sec. 


Current  analytical  models  for  estimating  the  local  rate  of  surface 
erosion  generally  assume  that  the  erosion  rate  Is  proportional  to  the  surface 
shear  that  would  exist  In  the  absence  of  lofting  (e.g..  Refs.  6,  7).  The  con¬ 
stant  of  proportionality  Is  deduced  from  Hartenbaum's  experimental  results. 

The  dependence  of  the  constant  of  proportionality  on  particle  and  fluid 
properties  Is  not  known.  Hence,  the  range  of  validity  of  this  approach  Is 
uncertain  with  departure  from  Hartenbaum's  test  conditions. 

The  present  study  Is  an  attempt  to  evaluate  local  erosion  rates  using 

conventional  turbulent  boundary-layer  theory.  The  case  considered  Is  one 

wherein  the  supply  of  surface  particles  Is  unlimited  and  the  particles  reach 

velocity  equilibration  with  the  ambient  fluid  soon  after  their  Injection  Into 

the  boundary  layer.  The  flow  Is  then  similar  to  the  case  of  a  turbulent 

a  q  1 A 

boundary  layer  with  surface  blowing.  *  *  It  Is  assumed  that  the  blowing 
rate  Is  such  that  the  surface  shear  is  reduced  to  the  threshold  value  of  shear 
at  which  surface-particle  mobility  Is  maintained.  A  similar  approximation  was 
made  for  the  saltation  regime  In  Ref.  2.  The  theory  Is  described  in  the  next 
section.  Turbulent-layer  properties  are  deduced  for  a  flow  of  air  over  a 
seml-lnf Inite  sand  bed  and  for  the  boundary  layer  behind  a  shock  moving  over  a 
sand  bed.  Comparisons  are  made  with  the  experimental  observations  of 
Hartenbaum^  and  Ausherman.^ 


II.  THEORY 


Integral  relations  for  a  turbulent  boundary  layer,  with  blowing,  are 
noted.  They  are  solved  subject  to  the  assumption  that  the  wall  shear  equals 
the  threshold  va^ue  at  which  surface-particle  mobility  Is  maintained. 
Expressions  are  derived  for  the  rate  of  dust  Ingestion  Into  a  turbulent 
boundary  on  a  flat  plate  and  behind  a  moving  shock  wave  (e.g..  Fig.  1). 

A.  BOUNDARY-LAYER  EQUATIONS 

Consider  a  turbulent  boundary  layer  with  surface  blowing  and  zero 
pressure  gradient.  The  effect  of  blowing  on  local  shear  and  boundary  layer 
thickness  can  be  approximated  by 


tn(l  +  B) 
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(la) 


and 
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Here  subscript  zero  denotes  nonblowing  values,  and  (p  v  )  denotes  the  local 

p  p  w 

rate  at  which  particles  are  Ingested  Into  the  boundary  layer.  Other  symbols 
are  defined  In  the  Appendix.  The  form  of  Eq.  (la)  can  be  deduced  from  a 

O 

consideration  of  Couette  flow,  with  blowing,  and  has  been  confirmed 

O 

experimentally  for  values  of  the  blowing  parameter  B  up  to  about  20. 

Equation  (lb)  is  deduced  from  the  momentum  Integral  equation  under  the 
assumption  that  the  ratio  6/6  is  relatively  independent  of  blowing.  Equation 
(lb)  Is  less  well-established  than  Eq.  (la). 


B.  THRESHOLD  SHEAR 

The  lowest  value  of  surface  shear  at  which  surface-particle  mobility  Is 
maintained  Is  termed  threshold  shear,  and  can  be  expressed  by 
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»  B  a  g  d  (2a) 

1  2 

where  B  is  a  dimensionless  constant  that  can  be  taken  equal  to  ' 

B  -  0.0064  (2b) 

2 

The  constant  B  may  be  viewed  as  the  ratio  of  shear  force  t  d  to  the 

3  ^ 

gravitational  force  a  gd  at  which  surface-particle  motion  can  be  maintained. 
C.  DUST-INGESTION  EQUATIONS 

We  now  assume  turbulent  flow  over  a  surface  consisting  of  an  unlimited 
number  of  small  particles.  We  also  assume  that  the  local  rate  of  ingestion  of 
particles  into  the  boundary  layer  is  Just  sufficient  to  reduce  the  local 
surface  shear  to  the  threshold  value.  Thus, 

Cf  -  '^f.c  = 


This  assumption  is  similar  to  one  made  in  Ref.  2  and  assumes,  in  effect,  that 
small  departures  of  from  ^  result  in  a  relatively  large  blowing  rate. 
Hence,  the  local  blowing  rate  corresponds  to  the  value  needed  to  reduce  the 
unblown  shear  C|  q  to  a  value  of  order  Also,  as  will  be  seen,  B  »  1 

for  cases  of  practical  interest.  For  these  cases  the  local  blowing  rate 
parameter  depends  on  the  log  of  C£.  Hence,  small  departures  of 

from  ^  should  not  significantly  affect  the  result,  and  Eq.  (3)  is 
applicable. 

If  the  surface  shear  in  the  absence  of  blowing  q  is  known,  the  local 
blowing  parameter  B  can  be  found  from  Eqs.  (la)  and  (3)  by  iteration,  namely 


The  local  blowing  rate  is  then 
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P  u 
e  e 


ln(l  +  B) 


(4a) 


(4b) 


and  is  seen  to  depend  only  weakly  on  B.  Other  boundary-layer  properties  are 
found  as  follows.  The  net  particle  flux  at  streamwlse  station  x  ■  L  can  be 
expressed  in  the  alternative  forms 


(5a) 


Equation  (3b)  assumes  velocity  equilibration  between  the  particles  and  the 
local  flow.  It  will  be  shown  that  for  the  present  applications,  Cj  q  ~ 

If  the  weak  dependence  of  £n(l  +  B)  on  x  is  ignored,  Eq.  (5a)  becomes 
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where  all  variables  are  evaluated  at  x  •  L.  In  order  to  evaluate  Eq.  (5b),  we 


.  1 _ c - 

u  ^Pn^u 

e  p  w 


That  is,  it  is  assumed  that  normalized  streamwise  velocity  and  particle- 
density  profiles  are  similar  and  have  a  power-law  variation  with  y/6 . 
Substitution  into  Eq.  (5b)  yields 


(p  )  u  6  ~  6  (I  +  n)(2  +  n) 
p  w  e 


(n  -  7)  (8b) 

Equation  (8a)  is  relatively  insensitive  to  n  and  the  value  for  n  ■  7  will  be 
used.  Equations  (lb),  (6)  and  (8)  yield  the  particle  density  at  the  wall. 


namely. 


^‘’p^w  45  *^f,0  B 

Pe  “  7  1  +  B 


and  the  particle  density  profile 
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The  particle  velocity  at  the  wall  is,  from  Eqs.  (lb),  (4b),  and  (9a) 


^'^p^w  7  6 


90  X 


(9c) 


These  equations  are  now  evaluated  for  the  case  of  incompressible  flow  over  a 
flat  plate  and  for  flow  behind  a  moving  shock  wave. 


D.  INCOMPRESSIBLE  FLOW  OVER  A  FLAT  PLATE 


Consider  Incompressible  flow  over  a  flat  plate  (Fig.  la).  In  the  absence 
of  blowing,  and  ignoring  surface  roughness  effects, 
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The  expressions  for  unblown  and  threshold  shear  become 
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Further,  consider  air  at  near-standard  temperature 

(Tg  - 

■  290  K). 
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Ing  numerical  values  then  apply: 

^e  ^e  -3 

(11a) 

195  7^  i.zia  X  lu 

3 

cm 

atm 

(^}  %  -  3.^15  X 

lo"^ 

cm 

sec 

(lib) 

e 

ogn  ^ 

II  a  1  R2q  X  10 

•4 

gr 

(lie) 

^  ^  J  *  I.OA7  X  L\J 

e 

sec  cm 

where  u  -  0.7.  For  sand-like  particles 
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where  and  x  are  In  atmospheres  and  centimeters  respectively  and  the 
dependence  on  Bd  Is  displayed.  Boundary-layer  properties  are  found  from 
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2.511  X  10 
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These  results  are  applicable  for  Reynolds  numbers  beyond  the  transition  value, 
Re^  >  10^.  Thus,  for  air,  the  region  of  applicability  is 


(M^Pe*)  >  ^ 


(13e) 


The  variation  of  boundary-layer  properties  with  Increase  In  q/C^  ^  Is 
Indicated  In  Table  1  and  Fig.  2.  For  the  large  values  of  Q/Cf  ^  of  present 
interest 


in(l  +  B)  = 
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Substituting  Eq.  (14)  Into  Eq.  (13)  provides  an  explicit  solution  for 

boundary- layer  variables  In  terms  of  the  parameter  o^^f,t*  ^  this  region, 

the  local  blowing-rate  parameter  (p  v  )  /p  u  and  the  boundary-layer  thickness 

p  p  w  e  e 

parameter  6/xhave  a  weak  (logarithmic)  dependence  on  q/C£  as  may  be 


Table 


Corresponding  Values  of  Blowing  Parameters  and  Shear  Ratio 

,o/Cf ,t 


fM 

B 

1.05  X  10® 

10" 

1.44  X  10® 

10® 

4.17  X  10® 

10* 

2.17  X  10* 

102 

1.45  X  10^ 

10^ 

1.09  X  10^ 

10* 

8.69  X  10^ 

io5 

7.24  X  10^ 

10® 

tn(l  +  B) 

9.53  X  io“2 
6.93  X  io“l 
2.40  X  10° 
4.62  X  10° 
6.91  X  10° 
9.21  X  10° 
1.15  X  10^ 
1.38  X  10* 


1  +  B 
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ln(l  +  B) 
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1.05  X  10° 
1.39  X  10° 
2.64  X  loO 
4.66  X  loO 
6.92  X  10® 
9.21  X  10° 
1.15  X  10* 
1.38  X  10* 


9.09  X  10“^ 

5.00  X  10“* 
9.09  X  10”* 
9.90  X  10“* 
9.99  X  10“* 

1.00  X  10® 
1.00  X  10® 
1.00  X  10® 


Variation  of  Turbulent  Boundary-Layer  Properties  with  Variation 
In  C*  n/C*  *  Cor  Seml-lnf Inlte  Plate  Plate  and  Moving-Shock 
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observed  in  Fig.  2.  This  result  is  consistent  with  the  previous  approximation 
that  the  shear  can  be  assumed  to  be  reduced,  by  blowing,  to  the  threshold 
value  for  surface  mobility. 

For  the  present  case  of  incompressible  flow,  -  p^.  The  dust  loading 
at  the  wall,  (p  )  /p  ,  is  approximately  equal  to  1.03  for  C,  n/Ce  ,  >  10.  The 

pWW 

wall  value  of  dust  loading  provides  an  upper  bound  on  the  dust  loading  within 
the  boundary  layer.  The  near -constant  value  of  (Pp)„/Py»  ^or  Cf ^q/Cj  ^  >  10, 
is  a  reflection  of  the  fact  that  the  blowing  rate  and  boundary-layer  thickness 
grow  at  the  same  rate  with  increase  in  q/C^ 

It  has  been  assumed  in  the  present  model  that  the  particles  remain  in 
suspension.  The  saltation  layer  study  by  Owen^  [e.g.  Eq.  (1)]  suggests  that 
this  assumption  may  be  invalid  when  q/C^  ^  ^  ^  (100),  B  <  0  (10^). 

However,  as  will  be  shown,  the  results  of  Hartenbaum  In  Table  2  and  Fig.  3a 
Include  data  at  o'^^f,t  “  ^1*7  which  are  consistent  with  the  present 
model.  Hence,  the  lower  limit  on  Cj  Q^^f  t*  which  the  present  results  are 

useful,  is  not  clear.  The  curves  In  Fig.  2  are  dashed,  for  <  50,  to 

reflect  the  uncertainty  in  the  lower  limit  of  Che  suspension  regime. 

E.  BOUNDARY  LAYER  BEHIND  A  MOVING  SHOCK 

The  boundary  layer  behind  a  moving  shock  is  illustrated  in  Fig.  lb.  A 
shock-fixed  coordinate  system  Is  used.  In  this  coordinate  system  the  wall  has 
a  velocity  Uy  equal  to  the  upstream  flow  velocity  Ug.  Flow  properties  up¬ 
stream  and  downstream  of  Che  shock  are  denoted  by  subscripts  1  and  e,  respec¬ 
tively.  Boundary-layer  parameters  are  normalized  by  the  free-stream  velocity 


relative  to  Che  wall,  Uy  -  u^.  Thus 

Cf,0  “  2  ^w,o/fPe('*w  ■  "e^^^ 

(15a) 

Cf.t  -  2  ^/IPe(“w  -  “e>^J 

(15b) 

B  -  2  P  P  w 

(15c) 

Cf  ^  P-(u  -  u  ) 
f,tew  e 
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The  normalized  velocity  and  particle-density  profiles  are 


u  -  u 
w 


u  -  u 
w  e 


1  - 


P  -  .  (1) 

p  V 


1/n 


(16) 


r.  •' 

^  'v*  ' 


The  net  particle  flux  at  x  ■  L  is  [assuming  ~  x“^^^  and  neglecting  the 

dependence  of  in(l  +  B)  on  x] 


5  *^f,0  B 


Pe(«„  -  “e>  ^  ® 


6o  1  +  B 


(17a) 


or,  alternately,  for  n  ■  7, 


. .  i 


M 

P  _  2W  +  7 

Tpn  (u  -  u  )  5  f2(W“-  1) 

p  w  we 


(17b) 


where  W  *  "  '*w^“e*  particle  density  and  velocity  at  the  wall  are, 

respectively. 


CJ 


_  5  72(W  -  1)  *^f.O 

“8  2W  +  7  1  +  B 


B 


and  [from  Eqs.  (lb),  (15c),  and  (18a)] 


^%^w  4  2W  -f  7  6 
“w  ■  "e"  ®  ^2(W  -  1)  7 


(18a) 


(18b) 


V  '.*■  \ 
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►  V  .• 

V*  Aj 


c>:y: 


In  the  absence  of  bloving,  Che  turbulent  boundary-layer  properties  can  be 
expressed  (for  air  with  T]^  ■  290  K  and  sand-like  particles) 


(Plx)^/5 


4.397  X  10-3  LW_t.(7/3)l 
(W  -  1)^'^ 


1/5 


xC 

~E 


fjO 

0 


0.0600 


W  +  (7/3) 
W  -  1 


S.O  _  1.991  X  10^  rw  -*•  (7/3)1^^^  (W  -  1)^^^ 

p,M  ^  ^f.t  6d/(6.4  X  10“^)  “ 

1  s 

The  solution,  with  dust  ingestion,  is  Chen  found  from 

^f,0  B 

Cj  t  ”  +  ®) 


(Pl*) 


1/5 


2.198  X  10 


_  ri!L: 

-3  '“W  + 


(W  -  1)^ 


1/5 


(p  V  ) 

P  P  w 

p  (u  -  u  ) 
e  w  e 


-  ln(l  +  B) 


^Pl^^  [W  -h  (7/3)1^^  §_ 

7.328  X  lo"^  (W  -  l)^'*^  * 


ln(l  +  B) 


1  2W  +  7  ^Pp^  w  _  B 
TITS’  W  +  (1/i)  “  1  +  B 

In  order  to  Insure  Chat  Che  boundary  layer  la  turbulent,  ve  require  the 

alenc  flat-plate  Reynolds  number^^  to  be  In  Che  range  Re  =  ~  1) 

X  e  e 

10  or,  for  air. 


(19a) 

(19b) 

(19c) 

(20a) 

(20b) 

(20c) 

(20d) 

equlv- 

■  */Pe  Z 


(21) 


(tt)  <M  p  X)  (W  -  1)^  >  4.4 

Equations  (20)  are  identical,  In  form,  with  Eqs.  (13)  and  are  included  in 

Fig.  2.  The  weak  logarithmic  dependence  of  the  blowing  rate  and  boundary- 
layer  thickness  parameters  on  Cj  g^^f  t  again  observed.  For  B  ♦  <•  and 
1  <  W  <  10  ,  the  particle  density  at  the  wall  is  in  the  narrow  range 
1.00  <  ^  1*23.  Evaluation  of  the  dust  loading  at  the  wall,  (Pp),,/Pj,* 

is  more  complex  than  for  the  previous  flat-plate  case.  If  we  assume  that  the 

wall  remains  at  its  initial  temperature,  >  Tp  ideal  gas  relations  indicate 


7  M, 


-  1 


(ideal  air) 


(22) 


1 


1 


The  dust  loading  at  the  wall,  for  ideal  air,  is  then 

W  +  (7/3) 


2.70 


6W 


2W  +  7 


7  Mf  -  1  '  "  ® 
s 


(23) 


For  strong  shocks,  (Pp)^7p^  falls  below  one.  This  is  due  to  the  fact  that, 

for  strong  shocks,  the  air  density  at  the  wall  rises  more  rapidly  with  shock 

strength  than  does  the  particle  density  at  Che  wall.  In  these  cases  the  peak 

dust  loading  (p  )/p  occurs  within  the  interior  of  Che  boundary  layer. 

P 

However,  this  feature  of  the  solution  may  be  an  artifact  of  the  simplified 
particle-density  profile  assumed  herein  [Eq.  (16)]. 


F.  VELOCITY  EQUILIBRATION 

It  has  been  assumed  that  Che  particles  Ingested  into  the  boundary  layer 
rapidly  achieve  velocity  equilibration  with  the  local  flow.  The  region  of 
validity  of  this  assumption  is  now  established. 

Assume  Chat  Che  dust  particles  leave  the  surface  vertically,  with  a 
velocity  V  The  initial  value  of  drag  coefficient  is,  for  spherical 


particles 


■  [(1/2)  (1,  (d/2)2] 


(24a) 


24 

Re. 


Rej  <  10 


(24b) 


0(1) 


(24c) 


where 


R«d 


(p.  V 


w  p,w 


Here  and  are  Quid  properties  evaluated  at  the  wall  [Eq.  (24b)]  corre¬ 
sponds  to  Stokes'  drag  law.  Let  t£  denote  the  time  for  the  particle  to  reach 
velocity  equilibration  with  the  local  fluid.  Let  denote  the  vertical 
distance  from  the  wall  at  which  velocity  equilibration  Is  reached. 

Approximate  values  for  these  quantities  are 


‘'E  "  y 

p,w  p 

where  (dVp/dt)„  ■  Fp/m.  For  spherical  particles 

X  ^  ^  g  X 

d  6  “  3  p  C_  ? 

w  D,p 


In  the  Stokes  region.  Re.  <  10, 

Q 


1 

2.02  X  10^ 


M  • 
e 


(plate) 


(25) 


(26) 


(27a) 


-  *  (shock) 


(27b) 


2.5  X  10 


1 


:  A 

The  quantity  ♦  equals  one  for  T,  -  290  K,  d  ■  10”^  cm,  and  0-2.5  gr/cm^. 

^  • 

For  Rej  >  10,  Eq.  (26)  is  evaluated  by  assuming  -  1.  In  this 

case,  -  Pg  for  a  flat  plate,  and  p^  is  evaluated  from  Eq.  (22)  for  a 
moving  shock. 

The  assumption  of  rapid  equilibration  is  valid  when 

«  1  (28) 

which  can  be  estimated  by  using  Eqs.  (26)  and  (27).  For  a  given  flow,  the 
assumption  of  equilibration  becomes  more  valid  as  x  (and  therefore  5) 

Increases. 

Equation  (25)  neglects  the  effect  of  cross  flow  on  This  neglect  is 

valid  in  the  Stokes  flow  regime  but  results  in  an  overestimate  of  t^/6  for 

£• 

Re  >  10  (Ref.  12).  The  latter  estimate  is  believed  adequate  for  present 
d 

purposes. 


III.  DISCUSSION 


Dust  Ingestion  Into  a  turbulent  boundary  layer  on  a  semi-infinite  flat 
plate  and  behind  a  moving  shock  has  been  studied  experimentally  in  Refs.  4  and 
5,  respectively.  These  experimental  results  are  compared  herein  with  the 
predictions  of  the  present  theory. 

A.  FLAT  PLATS 

Hartenbaum  has  reported^  boundary-layer  properties  and  dust-erosion  rates 
at  a  distance  of  x  *■  488  cm  downstream  of  the  leading  edge  of  a  silica  sand 
bed.  The  mean  particle  diameter  was  d  >  0.025  cm  and  the  test  gas  was  air  at 
standard  conditions.  Three  free-stream  velocities  (M^  ■  0.100,  0.212,  and 
0.336)  were  considered.  His  results  are  indicated  in  Fig.  3  and  Table  2.  The 
predictions  of  the  present  theory  are  Included  in  Table  2  and  Fig.  3a. 

The  local  erosion  rate  was  correlated  by  Hartenbaum  in  the  form  (Fig.  3a) 
(p  V  )  -  6.29  X  lo"^  u  (ibm/ft^  sec)  (29a) 

p  p  w  6 

where  u^  is  in  ft/sec.  For  the  range  of  flow  conditions  in  Ref.  4,  the 
present  theory  gives 

(P  “  2.70  X  lo’^  u  (Ibm/ft^  sec)  (29b) 

p  p  w  e 

The  present  theory  underestimates  the  local  erosion  rate  by  a  factor  of  0.4. 
The  exponential  dependence  on  u^  is  remarkably  close.  It  should  be  noted, 
however,  that  Eq.  (29b)  represents  a  tangent  to  the  corresponding  curve  in 
Fig.  2  and  therefore  is  applicable  only  for  a  limited  range  of  u^  and  x. 

Local  experimental  boundary-layer  thickness  can  be  estimated  from  the 
velocity  profiles  in  Fig.  3b  and  is  Included  in  Table  2.  The  predicted 
boundary-layer  thickness  is  larger  than  the  experimental  value  by  a  factor 
that  varies  from  2.0  to  1.3  as  increases  from  0.100  to  0.336.  The 
experimental  value  of  6,  however,  is  not  a  well-defined  quantity.  The  edge  of 


6  8 10^ 


6  8  10'" 


2  4 

Ug,  ft/sec 

Fig.  3.  Experimental  Data  at  x  -  488  cm  for  a  Turbulent  Air  Boundary 
Layer  over  a  Seml-lnf Inlte  Silica  Sand  Bed.  d  *  0.025  cm,  p, 
I  atm,  T_  -  290  K.  From  Kef.  4.  (a)  Local  mass-eroalon 


the  dust  layer  (see  Fig.  2c)  appears  to  extend  to  values  of  y  which  are 
greater  than  the  experinental  values  of  6  Inferred  from  the  velocity  profiles 
in  Fig.  2b  and  are  in  better  agreement  with  present  boundary-layer  thickness 
predictions.  The  experimental  velocity  profile  in  Fig.  3b  indicates  a  power 
law  u  ~  with  n  ■  1.6,  as  opposed  to  the  value  n  ■  7  used  herein.  The 

reduction  in  n,  with  an  increase  in  blowing  rate,  is  physically  realistic. 

The  present  theory  can  be  generalized  to  consider  velocity  profiles  of  the 
form  but  the  extension  does  not  appear  warranted  at  the  time. 

The  dust-density  profiles  presented  in  Ref.  4  have  been  extrapolated  to 
the  wall  in  order  to  provide  a  rough  estimate  of  (Pp)„/Pg*  Values  of  2  to  3 
are  inferred.  These  are  in  approximate  agreement  with  the  predicted  value  of 
1.03. 

The  theoretical  estimate  for  the  parameter  Z-/6  Increases  from  0.3  to  0.9 

£• 

as  Mg  increases  from  0.100  to  0.336.  Hence  the  applicability  of  the  theory  at 
the  larger  values  of  Mg  may  be  questioned.  However,  in  view  of  the  fact  that 
no  attempt  was  made  to  adjust  the  single  free  parameter  in  the  present  study, 
6,  the  agreement  between  theory  and  experiment  is  surprisingly  good.  The 
experimental  and  theoretical  predictions  of  local  erosion  rate,  boundary -layer 
thldcness,  and  dust  density  near  the  wall  agree  to  within  about  a  factor  of  2. 
No  systematic  divergence  between  theory  and  experiment  is  noted  as  1^/6 
increases  from  0.3  to  0.8.  This  si^gests  that  the  present  theory  may  be 
applicable  in  the  extended  range  i^/6  <  l.O,  as  defined  herein,  to  within  an 
accuracy  of  about  a  factor  of  2. 

B.  BOUNDARY  LAYER  BEHIND  A  MOVING  SHOCK 

Ausherman^  has  measured  the  turbulent  boundary  layer  behind  a  shodc 
moving  over  a  silica  dust  bed.  The  mean  particle  size  was  d  ■  0. 9  x  10~^ 
cm.  Measurements  were  made  at  a  fixed  station,  as  a  fmctlon  of  time  after 
shock  arrival,  for  shocks  moving  into  standard  air  at  Mg  ■  1.41,  1.74,  and 
2.25.  Time  after  shock  arrival  and  distance  behind  the  shock  are  related  by 
t  >  x/Ug.  The  variation  of  dust-cloud  height  with  time  was  evaluated  by 
observing  the  time  at  which  optical  extinction  commenced  at  each  of  four 
vertical  test  positions.  The  results  are  given  in  Fig.  4a.  The  variation  of 
dust  density  with  height  is  indicated  in  Figs.  4b  to  4d  as  a  function  of  time 
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after  shock  arrival  and  shock  Mach  number.  The  experimental  results  are 
compared  with  theoretical  predictions  in  Table  3. 

Comparison  of  the  theoretical  and  experimental  estimates  for  dust  cloud 
height  Indicates  that  the  theory  overpredicts  the  height  by  factors  that  range 
from  3.2  to  1.9  for  the  data  In  Table  3.  The  discrepancy  Is  reduced  with  an 
increase  In  x  and  M^.  The  magnitude  of  the  overpredlctlon  Is  about  the  same 
as  for  the  flat-plate  case. 

The  dust-density  profiles  were  extrapolated  to  provide  estimates  of  the 
dust  density  near  the  wall.  The  extrapolated  values  are  In  approximate 
agreement  (to  within  about  a  factor  of  2)  with  Che  theoretical  predictions. 

Values  of  Re^  and  are  Included  in  Table  3.  The  latter  has  been 

calculated  from  Eq.  (26)  by  using  Cq  p  >  1.  It  is  seen  that  ranges  in 

value  from  3.3  to  0.50,  %flth  the  lower  values  corresponding  to  Increased 

and  X.  These  values  of  i./d,  however,  are  believed  to  be  overestimates  that 

‘  12 
are  due  to  the  neglect  of  cross  flow  in  Che  non-Stokes  regime.  The 

agreement  between  theory  and  experiment  appears  again  to  be  within  about  a 

factor  of  2  for  those  cases  where  1-/6  <  1. 

E  ~ 
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IV.  CONCLUSION 


The  present  theory  contains  a  number  of  critical  assumptions.  The  most 
Important  are  that  (1)  conventional  turbulent  boundary-layer  correlations 
apply  for  very  large  values  of  the  blowing  parameter  B,  and  that  (2)  the 
effect  of  dust  Ingestion  Is  to  reduce  the  wall  shear  to  the  threshold  value 
needed  to  maintain  surface-particle  mobility.  In  addition,  it  has  been 
assumed  that  the  local  blowing  rate  Is  not  particle  limited  and  that  the 
particles  are  In  local  velocity-equilibration  with  the  ambient  fluid.  Simple 
normalized  profiles  were  assumed  In  order  to  calculate  dust  loading. 

The  theory  contains  a  single  free  numerical  parameter  6,  which  was  taken 

to  equal  the  value  recommended  by  Owen  in  his  study  of  the  saltation 
2 

process.  It  Is  somewhat  remarkable  that  the  predicted  erosion  rates  and 
boundary-layer  thicknesses  agreed  to  within  about  a  factor  of  2  with  the 
experimental  results  of  Hartenbaum  and  Ausherman.  This  agreement  may  be 
fortuitous,  since  the  amount  of  experimental  data  Is  meager  and  some  of  the 
data  (particularly  Ausherman' s)  violate  the  assumption  £  /6  «  I. 

Ci 

A  major  result  of  the  present  study  is  that  for  q/C^  ^  »  1,  the 
boundary-layer  parameters  have  a  logarithmic  dependence  In  C£  o^^f,t* 
hence  are  only  weakly  dependent  on  the  latter  parameter.  This  should  simplify 
the  estimate  of  local  erosion  rates  in  flows  of  Interest.  However,  further 
confirmation  of  the  present  theory  Is  needed. 
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APPENDIX:  SYMBOLS  AND  SUBSCRIPTS 


speed  of  sound,  cm/sec 

blowing  parameter,  Eq.  (1) 

particle  drag  coefficient,  Eq.  (24) 

local  shear  coefficient 

shear  coefficient  In  absence  of  blowing 

threshold  shear  coefficient,  Eq.  (3) 

particle  diameter,  cm 

drag  force  on  particle 

gravitation  constant,  cm/sec 

streamwlse  station,  cm 

particle-velocity  equilibration  distance,  cm 
Mach  number,  u/a 

net  particle  flux  at  station  L,  Eq.  (5) 
particle  mass,  gr 

velocity  profile  exponent,  Eq.  (7) 
number  of  particles  per  unit  volume 
pressure,  atm 
particle  Reynolds  number, 

Reynolds  number, 

equivalent  flat-plate  Reynolds  number,  (W 

temperature,  Kelvin 

fluid  velocity  In  x,y  directions 


x.y 

6 

e 


u 

p 


p 


p 


a 


T 

W 

ID 


particle  velocity  in  x,y  directions 

1  /2 

friction  velocity,  (2t^P^) 
density  ratio  across  shock,  P^/Pj 
streamv’lse  and  vertical  direction, 
boundary-layer  thickness 
boundary-layer  momentum  thickness 
viscosity,  gr/sec  cm 

fluid  density,  gr/cm^ 
particulate  density,  Opm 

single-particle  density,  gr/cvr 
surface  shear 

viscosity  lav  exponent,  Eq.  (11) 


Subscripts 

e  external  to  boundary  layer 

p  particle  value 

t  threshold  value,  Eq.  (2) 

w  surface  value 

0 


1 


nonblowing  value 
upstream  of  shock  wave 
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The  Aerospace  Corporation  functions  as  an  "architect-engineer**  for 
national  security  projects,  specializing  in  advanced  military  space  systems. 
Providing  research  support,  tne  corporation's  laboratory  Operations  conducts 
experimental  and  theoretical  investigations  that  locus  on  the  application  of 
scientific  and  technical  advances  to  such  systems.  Vital  to  the  success  of 
these  investigations  is  the  technical  staff's  wide-ranging  expertise  and  its 
ability  to  stay  current  with  new  developments.  This  expertise  is  enhtanced  by 
a  research  program  aimed  at  dealing  with  the  many  problems  associated  with 
rapidly  evolving  space  systems.  Contributing  their  capabilities  to  the 
research  effort  are  these  Individual  laboratories: 


Aerophysics  Laboratory:  Launch  vehicle  and  reentry  fluid  mechanics,  heat 
transfer  and  flight  dynamics;  chemical  and  electric  propulsion,  propellant 
chemistry,  chemical  dynamics,  environmental  chemistry,  trace  detection; 
spacecraft  structural  mechanics,  contamination,  thermal  and  structural 
control;  high  temperature  thermomechanics,  gas  kinetics  and  radiation;  cw  and 
pulsed  chemical  and  excimer  laser  development  Including  chemical  kinetics, 
spectroscopy,  optical  resonators,  beam  control,  atmospheric  propagation,  laser 
effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions, 
atmospheric  optics,  light  scattering,  state-specific  chemical  reactions  and 
radiative  signatures  of  missile  plumes,  sensor  out-of-f leld-of-view  rejection, 
applied  laser  spectroscopy,  laser  chemistry,  laser  optoelectronics,  solar  cell 
physics,  battery  electrochemistry,  space  vacuum  and  radiation  effects  on 
materials,  lubrication  and  sucfice  phenomena,  thermionic  emission,  photo¬ 
sensitive  materials  and  Infrared  detectors,  atomic  frequency  standards,  and 
environmental  chemistry. 

Computer  Science  Laboratory:  Program  verlf Icat Ion,  program  translation, 
performance-sensitive  system  design,  distributed  architectures  for  spacebome 
computers,  fault-tolerant  computer  systems,  artificial  Intelligence,  micro¬ 
electronics  applications,  communication  protocols,  and  computer  security. 

Electronics  Research  Laboratory:  Microelectronics,  solid-state  device 
physics,  compound  semiconductors,  radiation  hardening;  electro-optics,  quantum 
electronics,  solid-state  lasers,  optical  propagation  and  communications;  micro- 
wave  semiconductor  devices,  mlcrowave/milllmeter  wave  measurements,  diagnos¬ 
tics  and  radiometry,  microwave/millimeter  wave  thermonic  devices;  atomic  time 
and  frequency  standards;  antennas,  RF  systems,  electromagnetic  propagation 
phenomena,  space  communication  systems. 

Materials  Sciences  Laboratory;  Development  of  new  materials:  metals, 
alloys,  ceramics,  polymers  and  their  composites,  and  new  forms  of  carbon;  non¬ 
destructive  evaluation,  component  failure  analysis  and  reliability;  fracture 
mechanics  and  stress  corrosion;  analysis  and  evaluation  of  materials  at 
cryogenic  and  elevated  temperatures  as  well  as  In  space  and  enemy-induced 
environments « 

Space  Sciences  Laboratory;  Magnetospheric,  auroral  and  cosmic  ray 
physics,  wave-particle  interactions,  magnetospheric  plasma  waves;  atmospheric 
and  Ionospheric  physics,  density  and  composition  of  the  upper  atmosphere, 
remote  sensing  using  atmospheric  radiation;  solar  physics,  infrared  astronomy, 
Infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  ionosphere  and  magnetosphere; 
effects  of  electromagnetic  and  particulate  radiations  on  space  systems;  space 
instrumentation. 


